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3 ; ABSTRACT 

The 511 keV positron annihilation emission remains a mysterious component of the 
1^ I high energy emission of our Galaxy. Its study was one of the key scientific objective 

of the SPI spectrometer on-board the INTEGRAL satellite. In fact, a lot of observing 
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time has been dedicated to the Galactic disk with a particular emphasis on the cen- 



hH ' tral region. A crucial issue in such an analysis concerns the reduction technique used 



to treat this huge quantity of data, and more particularly the background modeling. 
0^. Our method, after validation through a variety of tests, is based on detector pattern 

Q . determination per ~ 6 month periods, together with a normalisation variable on a few 

hour timescale. The Galactic bulge is detected at a level of 70 a allowing more detailed 
C^ I investigations. The main result is that the bulge morphology can be modelled with 

two axisymmetric Gaussians of 3.2° and 11.8° FWHM and respective fluxes of 2.5 and 
5.4 xlO~^ photons cm~'^ s~^. We found a possible shift of the bulge centre towards 
^f-^ ■ negative longitude at I = —0.6 ± 0.2°. In addition to the bulge, a more extended struc- 

10 . ture is detected significantly with flux ranging from 1.7 to 2.9 xlO^'^ photons cm~^ s~^ 

^ . depending on its assumed geometry (pure disk or disk plus halo). The disk emission is 

r — [ , also found to be symmetric within the limits of the statistical errors. 
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Subject headings: Galaxy: general — Galaxy: structure — gamma rays: observations 



. F~ , 1. Introduction 

The detection of the 511 keV positron annihilation line emission from the central region of our 
Galaxy was one of the early and important successes of gamma-ray astronomy (Leventhal et al., 
1978). Although positron annihilation gives rise to the strongest 7-ray line signal from our Galaxy, 
three decades of dedicated observational and theoretical efforts failed to unveil the origin of the 
Galactic positrons. 



^Based on observations with INTEGRAL, an ESA project with instruments and science data centre funded by 
ESA member states (especially the PI countries: Denmark, France, Germany, Italy, Spain, and Switzerland), Czech 
Republic and Poland with participation of Russia and USA. 



With existing instrumentation, the emission appears to be diffuse; no point sources of annihila- 
tion radiation have yet been detected (Malet et al., 1995; Milne et al., 2000; Knodlseder et al., 2005; 
Teegarden et al., 2005; De Cesare et al., 2006;Weidenspointner et al., 2006). Its spatial distribution 
seems to be symmetric around the Galactic Centre (bulge), with an extent of 6 — 8° (Full Width at 
Half-Maximum) and a 511 keV flux measured to be around (1-3) xl0~^ photons cm~'^ s~^ (Lev- 
enthal et al., 1978; Milne et al., 2000; Kinzer et al., 2001; Knodlseder et al., 2005; Weidenspointner 
et al, 2008; Bouchet et al., 2008). 

There is no firm conclusion as to the origin of the positrons due to the modest angular resolution 
and limited sensitivity of 7-ray instruments. Further, the theoretical interpretation is complex since 
the distribution of the potential sources of positrons is unknown and many uncertainties remain 
on the composition of the interstellar gas, the structure and strength of the magnetic field in the 
Galaxy and the physics of positron diffusion and thermalization. 

Annihilation radiation from the disk and/or halo is still more difficult to study because of 
its lower surface brightness (Milne et al., 2000; Kinzer et al., 2001) but it potentially provides 
complementary clues to the positron production processes involved. Only a few instruments have 
been able to spatially resolve the disk emission from the brighter bulge one. Previous sparse 
measurements had shown that the annihilation emission from the disk is brighter in the longitude 
range [/| < 18 — 35° with very poor indications on the latitude extent (Gehrels et al., 1991; Milne 
et al., 2000; Kinzer et al., 2001). Using more than four years of data from the SPI/INTEGRAL 
imaging spectrometer, its spatial distribution has been constrained in the longitude range |/| < 100° 
and latitude extent |6| < 10° (Weidenspointner et al., 2008). 

In the following, we first present briefly the instrument and data selection. Then, we describe 
the analysis method developed for our purpose. The background determination being a critical 
issue for high energy data, particularly for extended emission and when the data encompasses a 
very long period, we deeply investigate this point and make extensive statistical tests to validate our 
results. Finally, the 511 keV sky distribution is derived using two approaches: an imaging method 
and a sky model fitting allowing to test several parameterized model distributions representing the 
Galactic bulge, disk and/or halo. All the results are discussed and summarized in the last part. 



2. Instrument and observations 

The ESA's INTEGRAL observatory was launched from Baikonour, Kazakhstan, on 2002 Oc- 
tober 17. The spectrometer SPI (Vedrenne et al. 2003) observes the sky in the 20 keV to 8 
MeV energy range with an energy resolution ranging from 2 to 8 keV. It consists of an array of 
19 high-purity germanium detectors operating around 80 K, the spectroscopic performance being 
maintained thanks to regular annealings of the detection plane. Its geometric surface is 508 cm^ 
with a thickness of 7 cm. The assembly is surrounded by a 5 cm thick bismuth germanate (BGO) 
which provides background rejection. In addition to its spectroscopic capability, SPI can image 



the sky with a spatial resolution of ~ 2.6° (FWHM) over a field of view of 30°, thanks to a coded 
mask located 1.7 m above the detector plane. The instrument's in-flight performance is described 
in Roques et al. (2003). Because of the small number of detectors, SPI's imaging capability relies 
on a specific observational strategy, which is based on a dithering procedure by steps of ~ 2°. 
Each pointing lasts between 30 and 60 minutes. We have analyzed observations recorded from 
2003, February 22 to 2009 January 2, covering the entire sky. Data polluted by solar flares and 
radiation belt entries are excluded. In order to completely suppress the effects due to radiation 
belts and to secure the data set, we suppress also systematically the first and last 10 exposures of 
each revolution. 

Finally, after image analysis and cleaning, ~ 1.1 x 10^ s corresponding to 39294 pointings (or 
exposures) are kept for the present study. The resulting exposure map is depicted in figure [TJ 



3. Data analysis 

The signal recorded by the SPI camera on the 19 Ge detectors is composed of contributions 
from each source (point-like or extended) in the field of view convolved by the instrument aperture, 
plus the background. For extended/diffuse sources, we assume a spatial morphology given by an 
analytical function or an emission map. For Ns sources located in the field of view, the data Dp 
obtained during an exposure (pointing) p, for a given energy band, can be expressed by the relation: 

Ns 

Dp = Y.Rp^oSp,3+Bp (1) 

i=i 

where Rpj is the response of the instrument for the source j, Spj is the flux of the source j, and 
Bp is the background recorded during the pointing p. Dp, Rp,j, and Bp are vectors of 19 elements. 
For a given set of Np exposures, we have to solve a system of Np equations (1). For that, it is 
mandatory to reduce the number of unknowns. 

A first way to do that takes advantage of the relative stability of the background pattern to 
rewrite the background term as: 

Bp = Ap X U X tp (2) 

where Ap is a normalization coefficient per pointing, U is the "uniformity map" or background 
count rate pattern on the SPI camera and tp the effective observation time for pointing p. U and 
t are vectors of 19 elements (one per detector). 

The number of unknowns (free parameters) in the set of Np equations is then Np x (Ng + 1) (for 
the Ng sources and the background intensities). Fortunately, a further reduction can be obtained 
since, at 511 keV, point sources are weak and constant within their error bars. Concerning the 



background intensity we consider it variable and test several timescales. Finally, we have Ns+Nf^ {Nf, 
being the number of time bins for the background) free parameters to be determined. Practically, 
Nb ^ 39000 if the background intensity is allowed to vary on the one exposure timescale (~ 2700 
s) and Nb ^ 6000 if the background intensity varies on a ~6 hour timescale. 

The global background count rate pattern on the SPI camera can also vary but on a much 
longer timescale (several months) . We did not find any noticeable pattern variation over 6 months of 
data, corresponding to the interval between 2 annealings. With this compromise: hours variability 
for the intensity parameter and modest (~ 6 months) variability for the background pattern , we 
obtain adequate fits to our measurements with an optimized number of parameters. 

We performed our analyses in the 508.25-513.75 keV band. Such a 5.5 keV wide band centered 
at 511 keV takes into account the Germanium energy resolution (FWHM ~ 2.05 keV ) including 
its degradation between two consecutive annealings (5%). At this energy, the gain calibration 
(performed orbit-wise) accuracy is better than ±0.01 keV. Finally, only single-detector event data 
are used (Roques et al., 2003). 



3.1. Core algorithm 

The tools used for background modeling, imaging, and model fitting were specifically developed 
for the analysis of SPI data and described in Bouchet et al. (2008). To determine the sky model 
parameters we adjust the data through a multi-component fitting algorithm, based on the maximum 
likelihood test statistics. We used Poissonian statistics to evaluate the quality of various sky models. 
The core algorithm to handle such a large, but sparse system is based on MUMPS g software 
together with an error bar computation technique dedicated and optimized for the INTEGRAL/SPI 
response matrix structure (Amestoy et al., 2006; Rouet, 2009; Tzvetomila, 2009). 



3.2. Background 

3.2.1. Modelling 

The modeling of the instrumental background is an important issue within the data analysis. 
Since 2 detectors failed (detector 2 on Dec. 7, 2003 and detector 17 on Jul. 17, 2004) during the 
period spanned by the observations, the distribution of the instrumental background in the detector 
plane changes significantly. 

This uniformity map or background pattern (eqn. 2) can be derived from empty field obser- 
vations. The dedicated INTEGRAL/SPI "empty field" observations are rare. But by making the 



^developed by the IRIT/ENSEEIHT laboratory http://mumps.enseeiht.fr 



hypothesis that the 511 keV emission is essentially concentrated along the Galactic plane and in the 
bulge, many exposures can be considered as empty-field observation at this energy. In practice we 
built a set of empty-fields using exposures whose pointing latitude satisfies \b\ > 30°. Furthermore, 
we exclude pointings closer than 30° from Cyg X-1 and the Crab Nebula since they are potential 
emitters in the 505-515 keV energy range. Finally, 5729 from the 39294 exposures can be used to fix 
the background patterns. Table 1 displays periods between consecutive annealings (~ 6 months), 
the revolutions spanned, the total number of exposures, the number of exposures used to build/test 
the "empty-field" and the statistics about the location in longitude and latitude of these exposures. 



3.2.2. Basic characteristics 

We define the quantity u to quantify the properties of the background pattern for each "empty- 
field" period, 

U 1 ^^" 

with Umoy = -TJ-JJ Yl ^(^) (3) 



^....^moy- ^,jf 



'moy 



d=l 



N^JJ is the number of working detectors (19,18 or 17), u and U are vectors of N'^-'^ elements, u 
represents the detector plane count rates normalized to the mean. 

Then, for each empty field period, we determine the background pattern as the time averaged 
vector u. Figure [2] displays the resulting pattern for two periods. The difference between the two 
patterns is essentially due to the failure of two detectors (2 and 17). 

To follow the pattern variation during the mission, we studied the evolution of this quantity 
for each detector. Figure [3] displays the evolution of the count rate measured on 2 detectors for 
each period using the empty field data set. This variation is very small, except at the detector 
failures. We also tried to improve the pattern determination by allowing a time dependency or by 
subdividing each "empty field" period in ~ 3 and -^ 1.5 month period but the F-test shows that 
the fit is not improved. In consequence, we used the ~ 6 month period interval as defined in table 
1 to determine the background pattern. 

As a final check, the distribution of residuals between our background model and the empty 
fields data has been tested against various statistical tests (Kolmogorov Smirnov, run test, moments 
of the distribution) and exhibit the expected normality properties. This distribution is plotted on 
fig. SI compared to a theoretical normal distribution. This confirms, a posteriori, the hypothesis of 
a constant pattern between successive annealings and the absence of any significant source. 

Another way to determine the background pattern is to solve eqn 1 for the source and the 
background simultaneously. While this method provides a better x^ (the global model being more 
flexible), the resulting patterns are model dependent. This is understandable in terms of an un- 
avoidable cross-talk between the extended source emission and the background pattern. As a 



consequence, the large scale emissions are partially suppressed and the statistical differences be- 
tween the models are reduced. Futhermore, the patterns obtained in such way are not perfectly 
compatible with the empty field data set. 



3.2.3. Background intensity variability 

The backgound model has been tested on our total dataset (39294 exposures, see table 1). The 
sky model consists of two point sources, namely Cyg X-1 and the Crab nebula, two axisymmetric 
Gaussians of FWHM 3.2° and 11.8° to describe the bulge and a Robin disk 1-2 Gyr (Robin et al., 
2003) to represent the spatial morphology of the disk. Eqn 1 is solved for Sj and Ap , U being 
fixed and computed from the empty-field dataset (5729 exposures) as described above. 

First, the background intensity Ap is assumed to vary on a timescale of ~ 6 hours. The residuals 
distribution between our model (sky plus background) and the data appears quite satisfactory 
(Fig. [5]). Secondly, Ap is allowed to vary on the pointing timescale (~ 1800-3600 s). The F-test 
shows that the second option does not improve the fit. In addition, the ~ 6 hour timescale analysis 
produces slightly reduced error bars (less free parameters) for background and source components 
while their intensities are perfectly compatible with those obtained with a more variable background. 

Finally, the standard configuration used for the subsequent analysis consists of a fixed bak- 
ground pattern per period between two annealings with an intensity variable on a 6 hours timescale. 



3.3. Investigation of the bulge morphology 

3.3.1. Maps 

The maps provide a model independent view of the sky but each pixel intensity is poorly 
determined due to the large number of unknowns. To gather reliable information, a 508.25-513.75 
keV sky map has been built using cells of size 61 = 5° x 5b = 5° (figures [6] and [7]) . 

Emission profiles are also derived using pixels of 61 = 16° x Sb = 16°, the central one being 
subdivided into a finer grid (pixel size of 51 = 3.2° x 5b = 3.2°). The reconstruction algorithm 
described in Bouchet et al. (2008) provides fiux and error for each sky pixel. In this energy 
range, the emission is almost entirely due to the annihilation line, the contribution of the other 
components being negligible. The emission profiles (figures [8] and [9]) are peaked towards and 
symmetric around the Galactic Centre. However, together with the map, they suggest that the 
emission is not limited to a bulge structure but exhibits an additional extended component (Galactic 
disk emission) revealed by the long exposure. 



3.3.2. Modelling the bulge with Gaussian profiles 

With more than 10^ s of exposure time, a single Gaussian (e.g. Knodlseder et al., 2005) does not 
provide an acceptable description of the bulge geometry. We obtain a better adjustment of the data 
with two Gaussians of FWHM 2.6 ° and 11.0 °, together with an extended component representing 
the Galactic disk modelled by a Robin 1-2 Gyr disk (~220° FWHM longitude extension, ~5.5° 
FWHM latitude extension). This two Gaussian model has been used for simplicity. Adding one 
or more sources simultaneously to the model, like GRS1758-258, GS1826-24 or H1743-22 known to 
emit above 100 keV, does not affect significantly the parameters of the Gaussians; none of these 
sources are detected above 2a. However, a significant change is obtained by introducing 1E1740.7- 
2942. It may reflect simply that the bulge profile contains small-scale diffuse emission structures, 
localized in the central ~ 3°. Furthermore the modest angular resolution of SPI makes it difficult 
to distinguish between point sources, point-like and small-scale diffuse emission. The Gaussian 
parameters depend weakly on the assumed disk geometry (for example the 240/i map is slightly 
peaked at the Galactic Center) but remain well inside the error bars. 



We also investigate a more complex bulge shape considering different widths in / and 6. This 
i A/ = (4.0lJ;^) and A6 = (2.5+J;°)° for the smallest Gaussian and M = (9.6+|^) and 



A6 = (15.2^3'q)° for the larger one. The statistics do not allow the derivation of any meaning- 
ful conclusion, but in all cases, these parameters stay compatible with those of the axisymmetric 
model. 



3.3.3. Possible shift of the Gaussian centroid 

In a following step, we let free the centroids of the Gaussians /qi ^o which are moving together. 
The fit gives ^o = (-0.64+°;^^)° and 5o = (0-06+[^;2o)°- These values depend slightly on the choosen 
configuration (in particular we tried several widths of the Gaussians), but all the parameters remain 
compatible. Fig.llOlshows the 1 a uncertainty zone in the longitude and the latitute for the centroid. 
In this configuration the introduction of an additional source is not needed anymore. 

To go further, we built a model composed of a set of nested shells of constant density centered 
at (/,6)=(-0.6,0)° with radii 0-3.0, 3-7, 7-10, 10-15 and 15-19°. We thus obtain the 511 keV flux 
distributions just assuming a radial emission. We split each shell into negative and positive part 
for the longitude and then for the latitude. Fig. [TT] and [12] display the radial profiles we obtained 
which appear quite compatible with the combination of the two off-centered axisymmetric Gaussians 
determined above. We thus conclude that the centroid emission appears slightly offset from the 
Galactic Centre direction. This result, while marginally significant, is all the more remarkable 
since Kinzer et al. (2001) have found a similar shift of the bulge towards negative longitude in the 
OSSE/CGRO data (1 a shift measured between -0.25 and -1°). 

Finally, we consider that the best configuration (hereafter called the reference model) con- 



sists of: 2 Gaussians centered at Iq = —0.6° and 60 = 0.0° representing the bulge plus an 
extended disk component (Robin 1-2 Gyr). Then the best fit FWHM of the Gaussians are 
(3.2+1;°)° and (11.8lJ|)° with fluxes of (2.48 ± 0.26) x 10"'' photons cm-'^ s'^ and (5.36 ± 0.5) x 
10~^ photons cm"'^ s~^ respectively, the disk flux being (1.83 it 0.27) x 10^'^ photons cm^"^ s^^, 
providing a Bulge to Disk ratio of 0.44. Using the ratio calculated by Knodlseder et al. (2005, Table 
3) to infer the 511 keV line luminosities and the factor (e"''/75ii) = 1.64 from Brown & Leventhal 
(1987), we obtain annihilation rates of 1.1 x lO^'^s"^ in the bulge and 0.8 x lO^^s^^ in the disk. 



3.4. Disk/halo emission 

We then attempted to describe the extended spatial distribution superimposed on the central 
bulge with various spatial geometries. Simple geometric shapes (i.e. two dimensional Gaussians) 
as well as more physical maps (CO, NIR, Robin disk (Robin et al., 2003)) were tested. 
Fig. [13] summarizes the result of the correlation map study, with the Ax^ ( which varies similarly 
as the reduced maximum log-likelihood ratio) plotted for each of the tracer maps. Best results 
are obtained for a Robin disk (in the 0.15-3 Gyr range corresponding to an old stellar population) 
or NIR/DIRBE 240/i and 1.25/x maps, which happen to be good tracers of the "^^Al line emission 
(Knodlseder et al., 1999). Simple bi-dimensional Gaussians of latitude FWHM ~ 5 — 7° and 
longitude FWHM ~ 250° give also good results. We note that the disk exhibits a larger longitude 
extension than the previoulsy reported values: first by OSSE (Kinzer et al., 2001), but the study is 
based on a longitudinally truncated data set distribution and more recently by Weindenspointner et 
al., 2008. Even though it is difficult to describe the emission in greater detail, this result represents 
a good indication for a bulge/disk structure. Extended disk structure flux (for the most plausible 
disks) is around 1.7 x 10~^ photons cnn~'^ s~^ and the Bulge/disk ratios range from 0.25 to 0.7. 

This analysis excludes single halo models (modelled by axisymmetric Gaussians). However, in 
the past, OSSE data have often been compared with bulge models including some non-Gaussian 
broadening (wings), featuring a bulge -|- halo central geometry. We have tested this kind of con- 
figuration by considering a stellar halo (or spheroid) model proposed by Robin et al. (2003). The 
central region profile is built following the law: 

N X {a/Ro)'^ for a > a^ (4) 

N X iajRorfora<ac (5) 

a2 = x^ + y2 + (z/ef (6) 

N represents a normalisation constant and x,y,z, the cartesian coordinates in the bulge reference 
frame. We obtain a good fit to the data with the axis ratio, e = 0.8, m = -2.6, ac=200 pc and 
Ro=8.5 kpc. The 240// and 1.25// NIR/DIRBE and Robin 1-3 Gyr maps remain the best tracers 
of the disk emission. We note further that assuming these geometries leads to higher B/D ratios 
(between 1.4 and 2.1 for the best models) since here B stands for the bulge plus the halo wings. 



3.4- 1- Disk symmetry 

To search for a disk asymmetry, we first used the galactocentric model and we divide the 
disk (Robin 1-2 Gyr) in 4 parts: a) -180° <1< -50°, b) -50° <1< 0°, c) 0° <1< 50° and d) 
50° <1< 180°. The fluxes of these different parts are respectively: (a) (0.84 ±0.27) (b) (0.44 ±0.07) 
(c) (0.29±0.07) and (d) (1.31 ±0.25) 10~^ photons C7n~'^ s~^. Those values, while compatible with 
those of Weidenspointner at al. (2008), support a symmetric disk within the quoted errors with 
a west to east ratio of 1.5±0.4. In fact, the main difference between both analysis concerns the 
background treatment. In Weidenspointner at al. (2008), the background pattern is determined 
once per orbit and simultaneously with the sky fluxes. It thus becomes model dependent, inducing 
additional uncertainties in the results (see discussion in section 3.2.2). 

Furthermore, we note that the disk flux in |/| < 50° depends on the assumed bulge geometry 
and position and is hence model-dependent. Now, using our best model for the bulge, i.e. the 
off-centered Gaussian described in 3.3.3, the split disk fluxes in region a) b) c) and d) are then 
(0.82 ± 0.27), (0.38 ± 0.07) ,(0.34 ± 0.07) and (1.32 ± 0.25) 10^^ photons cm'^ s'^ with a west to 
east ratio of 1.1 ±0.4. 

We have investigated the dependency of this ratio as a function of the bulge centroid. Fig. [H] 
displays the strong correlation between these 2 values and clearly proves that an asymmetric disk 
is not required by the data. 



4. Summary and conclusions 

We presented a detailed investigation of the morphology of the Galactic emission in the narrow 
annihilation line at 511 keV. The period covered by the data encompassing more than 6 years, the 
background determination is a key point in the analysis. We have chosen to describe the background 
with a series of fixed patterns (about one per 6 months), together with a normalisation adjusted by 
timebins of 6 hours. The patterns have been measured in the sky regions outside the Galactic disk, 
with the advantage to be independent of the sky model. This method has been validated through 
various statistical tests and represents a major difference between our analysis and that presented 
in Weidenspointner (2008). 

While the reality might be more complex, the bulge morphology is well described by two 
axisymmetric Gaussians, centered at Iq = —0.64° and bo = 0.06°, with FWHM of 3.2° and 11.8° 
and fluxes of (2.48±0.26) and (5.36±0.50) x 10~^ photons cm~'^ s~^ respectively. While marginally 
significant, it is similar to the shift of the bulge towards negative longitudes seen in the OSSE data 
(Kinzer et al., 2001). 

We confirm also the existence of a disk emission with a longitude extent over 150 — 250° and 
latitude over 5 — 7°. 
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Many spatial distributions can be used to model the disk/halo geometry. Among them, the 
best fits are obtained with disk distributions extending over ~ 200° x 5 — 6° (gaussian shapes, 
240;U and 1.25/Li NIR/DIRBE maps or Robin 1-3 Gyr models). The disk flux is around 1.7 
xlO~^ photons cm~'^ s~^ leading to a bulge-to-disk flux ratio of 0.4. This one ranges from 0.25 
(for the 1.25/i NIR/DIRBE map which presents a broad bump emission in the Galactic Centre 
region) to 0.7 (for less extended distributions like 150° x 5° gaussian). This kind of disk emission 
morphology, widened in longitude, suggests an old stellar population as the main Galactic positron 
source. 

When allowing a halo contribution (Robin stellar halo model), the flux reaches a value of ~ 
2.9 X 10~^ photons cm,~'^ s~^, with bulge-to- disk flux ratios evolving from 1.7 to 2.1 with the 
same trend than above. 

With our best flt determination, the disk is found to be symmetric with a west to east ratio 
of 1.1. It is clear that the disk results are very sensitive to the assumed bulge morphology, with a 
strong coupling between the parameters of both structures. Fig. [TH illustrates that, considering the 
error bars on the bulge position, it is impossible to claim a disk asymmetry. The observed trend 
can be understood in terms of unavoidable cross-talk between both emissions (bulge/disk) in the 
common region: when the bulge is allowed to move, the reconstructed disk flux is affected. 

The INTEGRAL mission has significantly advanced positron astronomy. We begin to have a 
reliable view of the annihilation emission from the inner Galaxy, and can hope to possibly constrain 
the extended emission. It should bring some breakthrough in the quest for the mysterious origin of 
the positrons. At the same time, spectroscopy of the annihilation emission has yielded important 
insights into the conditions of the medium in which the positrons annihilate (Churazov et al., 2005, 
Jean et al., 2006). Our work can be summarized as follows: 

• 511 keV emission is significantly detected ( 70o" in data space) towards the galactic bulge 
region, and, at a lower level (25cr) from the Galactic disk 

• There is no evidence for a point-like source in addition to the diffuse emission, down to a 
typical fiux limit of 4 x 10~^ photons cm'^'^ s~^ (^c). 

• Among the variety of possible disk geometries, we favor extended (~ 200° x 5 — 6°) disk 
distributions, suggesting an old stellar population as the main Galactic positron source. The 
bulge-to-disk flux ratio ranges from 0.25 to 0.7. 

When allowing a halo contribution, the bulge-to- disk flux ratios evolves from 1.7 to 2.1. 

• For our reference model (2 gaussians bulge + Robin 1-2 Gyr disk), we obtain annihilation 
rates of 1.1 x lO^^s^^ in the bulge and 0.8 x lO^^s"^ in the disk. 

• The disk asymmetry reported in Weidenspointner et al. (2008) is not supported by our 
analysis. 
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The bulge emission appears almost spherically symmetric around the Galactic Centre with 
an extension of 12°. Its centre seems slightly shifted towards negative longitudes at / = 
-0.64 ±0.20°. 
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Fig. 1. — 508.25-513.75 keV INTEGRAL SPI exposure map. Units are in cm^ x s. This map takes 
into account the differential sensitivity of SPI across its field of view. 
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Fig. 2. — 508.25-513.75 keV relative pattern for the empty-field corresponding to revolutions 97-131 
(dash-dot line) and 647-713 (solid line). 
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Fig. 3. — 508.25-513.75 keV pattern for detector 3 (dash-dot line) and 12 (solid line) as a function 
of the empty-field number. The dotted lines indicate when detector 2 and 17 failures occured. 
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Fig. 4. — 508.25-513.75 keV residual distribution for the entire empty field dataset (residuals 
minimum and maximum are respectively -4.02 and 4.46). The distribution is compared to a normal 
distribution with a variance of 0.94. The Kolmogorov-Smirnov test gives a distance d=0.005 and a 
probability p=0.94. The run test probability is 0.38 
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Fig. 5. — 508.25-513.75 keV residual distribution for our reference sky model (residual minimum 
and maximum are respectively -4.97 and 4.36. The distribution is compared to a normal distribu- 
tion having a variance of 0.996. The Kolmogorov- Smirnov test gives a distance d=0.0007 and a 
probability p=0.999, the run test probability is 0.05. 



17 





1 






1 i 










1 






1 


5 ' 










J 














■ 


■^H^^. 






B ■ 




I 




1 






■ 


^^^^^ 






t ■ 




I 


^1 


i 




■ 












■ 








u 




^H 


""\ 




















1 
















I^H 












1 


■ 








1 






^^^^^H 


^V 




■ 






■ 


■w^„. ._ 












r 




i 


1- ^^^^1 








p 


4^ 


^H 


^^^^ 


^^^^^^ 




r— 1 


-^H 1 






H T|| 




f 


^'^ 




••r^t 


■g 




la-i |-| 


^ 




60 jO 


i-i 


^^ 


n 


■ 


T-"-^ 




' 




*J , —r-\ 


■UM 


■^ 


rl 


-^ 


■ 


""-^"- 


-^ — ~^ 




■ ■ ^ 




=F=^ 1 1 


-J— . , 


3 




1 
1 


m ' 














^ 


1^^ 




J 




m 


1 

■ 


*. 














- 




1 


■^ 


I 


■■ 


W 














B i 






^^H 


_ 






i. 














E J 






^^1 




1 


^m 




1 


1 






r- ■ 


' 






T 


F^^ 


■ 


1 












■ 




■ 








■ ■ 
















3^ 


□ 














1 




trri 


H 






















i 











GcLxilio longhude 



Fig. 6.— 508.25-513.75 keV INTEGRAL SPI smoothed (top hat of 2 pixels) significance map. 
Pixel size is 5° x 5°. 
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Fig. 7.— 508.25-513.75 keV INTEGRAL SPI smoothed (top hat of 2 pixels) intensity map in 
photons cm^'^ s^^. Pixel size is 5° x 5°. 
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Fig. 8. — 508.25-513.75 keV longitude profile for \b\ < 8°. Dotted-dashed line corresponds to the 
sum of 3.2° and 11.8° axisymmetric Gaussians. Dotted line corresponds to the extended distribution 
(Robin 1-2 Gyr disk model). The sum of both has been integrated on the same bins as the data 
(dashed histogram). The sky is divided into 16° x 16° pixels except the central one which is divided 
into pixels of 3.2° x 3.2° 
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Fig. 9. — Same as fig. 8 but latitude profile for |/| < 8°. 
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Fig. 10. — 2D la error box on the bulge position for 3 disk models. The bulge is modelled with 
3.2° and 11.8° axisymmetric Gaussians. The disk is modelled by (a) a Robin 1-2 Gyr disk (solid 
line), (b) a Robin 1-2 Gyr disk segmented in 4 parts (short-dashed line) and c) a Gaussian with 
a longitude FWHM of 250° and latitude FWHM of 6° (long-dashed line) The best centre is found 
at (a) Iq = -0.64 ± 0.20 and 6o = 0.06 ± 0.20, (b) Iq = -0.57 ± 0.24 and 6o = 0.06 ± 0.20 and (c) 
Iq = -0.63 ± 0.20 and 6o = 0.06 ± 0.20. 
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Fig. 11.— 508.25-513.75 keV radial profile for b < 0° (filled circles) and 6 > 0° (diamonds). 
Dotted-dashed line corresponds to the sum of a 3.2° and 11.8° axisymmetric Gaussians. 
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Fig. 12.— 508.25-513.75 keV radial profile for I < -0.6° (filled circles) and / > -0.6° (diamonds) 
Dotted-dashed line corresponds to the sum of a 3.2° and 11.8° axisymmetric Gaussians. 
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Fig. 13. — Maximum log likelihood-ratio (MLR) as a function of the tracer map. The left part 
corresponds to physical maps (^ CO map (Dame et al., 2001)). The central part corresponds to 
pure halo emission modelled with axisymmetric Gaussians. In the last part, the disk is assumed 
to have a Gaussian morphology whose latitude and longitude FWHM are indicated horizontally 
and vertically respectively. Gaussians with longitude > 250° FWHM are not shown as they do not 
improve the MLR. The bulge is always modelled with the two off-centered Gaussians. 
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Fig. 14. — Top panel: MLR value as a function of the bulge centre position in longitude. Bottom: 
ratio of flux in— 50° <1< 0° to 0° <1< 50° as a function of the bulge centre position. The curves 
depict the results for two models: a) The bulge is modelled with two axisymmetric Gaussians of 
3.2° and 11.8°, the disk with a 1-2 Gyr Robin disk (solid line and filled circle), b) The longitude and 
latitude FWHM of the Gaussians modeling the bulge have been fitted for each longitude position. 
The disk is modelled by a bi-dimensional Gaussian of FWHM 250° x 6° (Dotted line and diamonds) 
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Table 1. Basic empty- field dataset construction. 



Inter annealing 


Empty-field limits 


Data total 


Percentage of 


Empty-field exposures distribution 


period number 


First-last revolutions 


exposure number 


empty-field 


/ < 071> 0°/b< 0°/b> 0° 


1 


44- 


- 92 


3328 


12.6 


192 / 226 / 8 / 410 


2 


97- 


131 


1816 


8.4 


58 / 95 / 101 / 52 


3" 


137- 


- 139 










140- 


- 204 


4334 


5.2 


128 / 97 / 129 / 96 


4" 


210- 


- 214 


277 


8.7 


/ 24 / 24 / 




215 ■ 


- 277 


3730 


13.0 


218 / 268 / 181 / 305 


5 


283- 


- 325 


2767 


18.8 


199 / 321 / 72 / 448 


6 


331 ■ 


- 394 


3550 


10.6 


138 / 240 / 186 / 192 


7 


401 ■ 


- 445 


2668 


15.2 


296 / 110 / / 406 


8 


452 - 


- 505 


3093 


15.2 


278 / 191 / 401 / 68 


9 


512 - 


- 564 


3309 


8.0 


51 / 213 / 46 / 218 


10 


571 ■ 


- 640 


4364 


19.2 


342 / 495 / 261 / 576 


11 


647- 


- 713 


3361 


25.0 


425 / 414 / 252 / 587 


12 


720 ■ 


- 775 


2697 


36.9 


877 / 117 / 407 / 587 



Note. — " detector 2 failure on Dec. 7, 2003 (between revolution 139 and 140) , ' detector 17 failure on Jul. 17, 
2004 (between revolution 214 and 215). There is enough empty-field exposures to compute the detector pattern before 
(revolutions 210-214) and after (revolutions 215-277) the failure of detector 17. In the case of detector 2 failure, the 
detector pattern can be derived only after the revolution 140 (revolutions 140-204). 



